Rubidium 87 Bose-Einstein condensate in an optically plugged quadrupole trap 
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We describe an experiment to produce *^Rb Bose-Einstein condensates in an optically plugged 
magnetic quadrupole trap, using a blue-detuned laser. Due to the large detuning of the plug laser 
with respect to the atomic transition, the evaporation has to be carefully optimized in order to 
efficiently overcome the Majorana losses. We provide a complete theoretical and experimental study 
of the trapping potential at low temperatures and show that this simple model describes well our 
data. In particular we demonstrate methods to reliably measure the trap oscillation frequencies and 
the bottom frequency, based on periodic excitation of the trapping potential and on radio-frequency 
spectroscopy, respectively. We show that this hybrid trap can be operated in a well controlled regime 
that allows a reliable production of degenerate gases. 
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I. INTRODUCTION 

Atom traps combining magnetic and optical potentials 
have been used since the achievement of the first Bose- 
Einstein condensates (BECs) in 1995 They take ad- 
vantage of the large trapping volume offered by mag- 
netic quadrupole traps which facilitate the loading of a 
large atom number from magneto-optical traps (MOTs). 
Moreover, they allow for efficient evaporative cooling dy- 
namics thanks to the initial linear shape of the con- 
finement, while still offering large trapping frequencies 
at the end of the evaporation. Furthermore, they pre- 
vent the trapped atoms from undergoing Majorana spin 
flips 0, relying on the dipole force induced by a laser 
beam to push the atoms outside the low-magnetic-field 
region where these losses are the highest. 

While solutions relying on red-detuned lasers have 
been demonstrated (3], an alternative consists in using 
a blue-detuned laser as an optical plug at the center of 
the magnetic quadrupole trap; see Fig. [1] The latter 
strategy has the advantage of minimizing resonant scat- 
tering processes responsible for atom losses and heating. 
A widespread choice for the design of the optical plug 
consists in using standard lasers at 532 nm. These lasers 
indeed commonly reach several watts of power and have 
long demonstrated their stability and reliability. This de- 
sign is particularly well suited for ^"^Na (sodium) atoms, 
given the vicinity of the D lines at 589 nm. It has hence 
recently allowed for the fast preparation of large ^'^Na 
BECs 11, i]. 

Here we show that the same setup can also be used 
to efficiently obtain ^^Rb (rubidium) BECs despite the 
large detuning of the optical plug laser from the ^^Rb 
main transition at 780 nm. Indeed, we are able to keep at 
all times a high collision rate while minimizing Majorana 
losses by modifying the trap shape during the evaporative 
cooling ramp. We typically obtain 2 x 10^ atoms in a 



quasi-pure BEC after a total evaporation time of 20 s. We 
finally present an experimental investigation of the trap 
parameters and describe how to characterize the cold gas 
during the cooling process. 

The paper is organized as follows: In Sec. |TT] we per- 
form a thorough analysis of the optically plugged trap. In 
Sec, mil we show how to optimize the evaporative cooling 
dynamics. In Sec. IIVI we give an overview of the experi- 
mental methods allowing us to infer the trap parameters 
and also explain how to characterize the cold gas along 
the way to a BEC. Additional details on the experimental 
setup are given in the Appendix. 
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Figure 1: (Color online) Principle of the optically plugged 
quadrupole trap. A blue-detuned focused laser beam plugs 
the central part of a quadrupole magnetic trap, where the field 
is zero. Left: The magnetic field is produced by a pair of coils 
of axis z aligned with gravity, with opposite currents. The 
beam propagates along the y axis. Right: The isomagnetic 
lines are represented by dashed lines in the xz plane, together 
with a line of constant intensity of the plug beam, orthogonal 
to the plane. 
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II. OPTICALLY PLUGGED TRAP 

In this section, we present a detailed description of 
optically plugged quadrupole traps, with the geometry 
depicted in Fig. [1] When the plug beam is centered ex- 
actly at the magnetic field zero, we are able to derive 
analytical expressions for the trap position and its cor- 
responding oscillation frequencies, as well as to quantify 
the sensitivity of these quantities to the relevant experi- 
mental parameters. These results are generalized to the 
case of a slightly off-centered plug beam. This study al- 
lows us to infer the dependence of the trap characteristics 
on the experimental parameters and in turn to estimate 
the heating rate due to their fluctuations, in particular 
regarding beam pointing stability. 



A. Trap geometry 
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As stated in the Introduction, optically plugged 
quadrupole traps result from the combination of a 
quadrupole magnetic field with a blue detuned laser beam 
focused near the symmetry center of the quadrupole, 
where the magnetic field is zero. Taking into account 
the gravitational field as well, the corresponding poten- 
tial can be written 



U{Y) = UB{v) + UD(v) + Mgz 



(1) 



where Ub(j) is the magnetic component of the potential, 
Ud{y) its optical component. A/ the atomic mass, and g 
the gravity acceleration. 

Depending on the relative choice of the quadrupole 
magnetic field's orientation to the plug beam direction, 
different trap configurations are possible: there can be 
a single minimum Q , a pair of minima [H 0] or an in- 
finite number of degenerate minima spread over a cir- 
cle Q. In the following we will consider only the case of 
a plug beam propagating along one of the horizontal axes 
of a magnetic quadrupole field of vertical symmetry axis, 
which corresponds to our experimental configuration (see 
Fig. [1]) . In general such a potential presents two minima 
separated by two saddle points; see Fig. O 

As depicted in Fig. [1] the quadrupole magnetic field 
is produced in our experimental setup by a pair of coils 
with vertical axis z and can be written 



-4z2 



(2) 



where a = fib' /h and fi ~ mpgpfis, with b' the horizon- 
tal magnetic gradient, gp the Lande factor in the ground 
state F, fiB the Bohr magneton, and mp the atomic spin 
projection on the magnetic field axis. The potential U b 
is hence cylindrically symmetric around z. 

The plug beam is produced by a green laser at 532 nm 
propagating along the y direction and focused to a waist 
Wq near the zero-field position. Denoting by {Xc,Zc) the 
focus position of the beam and neglecting the effect of 



Figure 2: Calculated trapping potential in the xz plane in the 
F — l,mF = —1 ground state resulting from the combination 
of the quadrupole field, the optical plug centered at position 
{xc,Zc) and gravity, see Eq.([T]). The parameters are given in 
Table IT] The lines are isopotentials separated by 20 kHz. The 
bold lines refer to the isopotential 1 kHz above the right min- 
imum at {xq,zq). The eigenaxes of the trap {x',z'), making 
an angle 6 with the reference axes, are also shown. Inset: In 
situ absorption image of ultracold atoms confined in the two 
minima of the optically plugged trap. 



the finite Rayleigh length — typically much larger than 
the atomic sample — the dipolar potential reads @ 



Ud (r) = Uo exp 



(X - Xc)'^ + iz~ Zc)'^ 



(3) 



where Uq is the maximum light shift at the beam waist. 
L^o is proportional to the plug power P and to the inverse 
squared waist, Uq esc P/wq [7[. 



B. Trap characteristics 

From Eqs. © and we see that the trap potential 
U is completely determined by five parameters, namely, 
6', P, wq, and {xc,Zc)- In this section, we restrict our 
study to the case of a centered beam, Xc = Zc = 0. Nev- 
ertheless, our main conclusions remain valid if the plug 
beam is slightly off centered, the trap positions and the 
corresponding oscillation frequencies being only weakly 
affected. 

Symmetries of the potential U allow us to make a few 
general statements about the trap geometry. First of 
all, the dependence oi U on y is restricted to the term 
Ub, even in y, which imposes the requirement that all 
the potential minima (and saddle points) should belong 
to the y = plane. Hence, the y axis will always be 
an cigenaxis of the trap in a harmonic approximation 
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around one of the trap minima. Second, U is even in x. 
As a consequence, the two saddle points belong to the 
X = plane, and the two trap minima are symmetric 
with respect to this plane. In the following we will hence 
concentrate on the minimum located on the right side 
(xo,0, zq). Because of gravity, the trap depth is limited 
by the lower saddle point belonging to the z axis, at 
position (0,0, Zs), with Zg < 0. 

Minimizing Ec^. Q , we can express the coordinates xq 
and zq as 



2^0 



roVl — 4e and 



(4) 



Here, we have introduced the dimensionless parameter e 
describing the relative strength of gravity and the mag- 
netic gradient as £ = [M g / {3ha)]^ . The effective radius 
To ~ -^/xq + 4zq can be written 



ro 



Wo, 



(5) 



where the dimensionless parameter ^ is the solution of 
the equation ^exp[— 2^^] = hawoV^ — 3e/[4C/o], verify- 
ing ^ > 0.5 Q. The isomagnetic surface Bq on which the 
minima of potential U are located is directly linked to rg 
through the relation 



(6) 



From Eq. ^ and Q, we see that a few conditions 
must be fulfilled to ensure the existence of a minimum 
at position (xo,0, zq): namely, e < 1/4 and Uo/[hawQ] > 
— where e = exp[l]. The first relation is usu- 
ally verified since with typical experimental parameters 
gravity is small compared to the magnetic field gradient. 
The second one requires that the light shift overcome the 
effect of the magnetic gradient on the size of the waist. 

We now perform a harmonic approximation of the po- 
tential at the position of the minimum (a;o,0,zo) a-nd 
derive the corresponding oscillation frequencies: 



ha 



la;,, 



(7a) 

(7b) 
(7c) 



The frequency iOy appears as the natural frequency scale 
of the trap and is generally the smallest one. Interest- 
ingly, from Eqs. © and ([Ta]) . we see that uiy is completely 
determined by the value of the magnetic field at the trap 
minimum, Bq, and the magnetic field gradient b' through 
the relation ujy = y/ jib'^ / M B^ . 

The two largest frequencies ujx' and w^/ correspond 
to the eigenaxes x' and z' which are rotated with re- 
spect to th e X and z axes by the angle 9, defined by 
tan^ = —\/e/ (1 — 4e) = Zq/xq (see also Fig. [2]). The an- 
gle 9 is usually small. We point out that the eigenaxis x' 
coincides with the line linking the center of the magnetic 
field (0,0) to the center of the trap {xo,zq). In a basis 



Parameter 


U 


P 


Wo 




Zc 


Uo/kB 


depth 


Unit 


G-cm-i 


W 


fj,m 






fiK 


fiK 


Value 


55.4 


5.8 


46 


5.5 


-0.5 


100 


5.3 


Uncertainty 


±0.6 


±0.1 


±3 


±3 


±3 


±15 






e 




C 


ro 




Uy 




Unit 








fim 


Hz 


Hz 


Hz 


Value 


0.0337 


1.47 


0.131 


77 


220 


76.6 


121.1 


Uncertainty ±0.0007 ±0.04 ±0.007 


±1 


±10 


±0.4 


±0.5 



Table I: Set of parameters used in the experiment, and char- 
acterized in Sec. IIVI leading to the isopotential lines plotted 
in Fig. [21 Atoms are prepared in the _F' = l,mF = — 1 ground 
state. &', P, and the oscillation frequencies Vi are measured 
(see Sec. IIV B)) : the other parameters and the related micer- 
tainty are deduced from these measurements. Here, Vi stands 
for uji/{2n), where i = x',y,z'. 



centered in (0, 0) and rotated by 6, the coordinates of the 
center of the trap are x'q = a/xq + Zq = ^/l ^ 3e tq and 
4=0. 

The positions of the saddle points as well as the trap 
depth can also be calculated analytically. The exact ex- 
pression is somewhat complicated. However, we notice 
that the distance of the saddle point from the center is 
of the same order as the distance xq or rp. The contri- 
bution of Ud is then about the same at both positions, 
whereas the contribution of Ub is twice as large at the 
saddle point, due to the larger gradient along z. Includ- 
ing the effect of gravity, we can thus estimate the trap 
depth to be haro — Mgro = haro(l — S^/e). 



C. Sensitivity to the trap parameters 



Magnetic gradient, beam power and waist 



From Eq. © and setting C = 1/(4^^ - 1) = ^y/^l', 
we can express the dependency of the effective radius vq 
on the experimental parameters as 



ro (X 



b'- 



ro (X 



1-3C 



(8) 



With our experimental parameters, presented in Table U 
the three exponents arc —0.26, 0.13 and 0.6 respectively. 
This implies that the effective radius rg is relatively stable 
to fluctuations of the plug beam power and the magnetic 
field gradient and is only slightly sensitive to the fluctu- 
ations of the beam waist, the latter being certainly one 
of the most stable experimental parameters. 

Fluctuations in ro are responsible for fluctuations of 
the trap position (xo,zo), which leads to linear heating 
through dipolar excitation, with a constant temperature 
time derivative T @. T is proportional to the power 
spectrum density (PSD) of the position noise, which can 
be linked to the PSD of the noise in b', P and wq at the 
trap frequencies through Eq. ([S]). Similarly, fluctuations 
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in ro also lead to fluctuations in the oscillation frequen- 
cies, which produce an exponential heating due to para- 
metric excitations, with a rate rparam proportional to the 
PSD at twice the trap frequencies [9[. From the relations 
([8]) , we can estimate the maximum allowed density power 
spectrum of relative fluctuations of 5', P and wo that en- 
sure a linear heating rate T and an exponential heating 



rate Fn 



below an arbitrary given threshold. 



As the trap minimum always lies in the plane y ~ Q, 
fluctuations in the parameters do not change this coor- 
dinate and the contribution of the y direction to dipolar 
heating is zero [l^l- On the other hand, as Zq = for a 
centered plug beam, the relative fluctuations of P and wq 
do not induce dipolar heating along z' . Fluctuations in b' 
do change z'q because of a change in the angle 6', through 
5zq = \J x\-\- z^ S9. However, even for b' , the dominant 
contribution to dipolar heating is along the direction x' . 

For a threshold at T = lnK-s~^, wc find that the 
maximum power spectral density of relative fluctua- 
tions should be 5'&'(i'x') < -100dB-Hz"\ Sp{vx') < 
-95dB-Hz-i, and S^„{i^:^>) < -llOdB-Hz-i for 6', P, 
and Wo, respectively. These requirements are easily ful- 
filled with commercial power supplies and lasers; see the 
Appendix. For parametric heating, all three directions 
have to be taken into account. If the noise spectrum 
at frequencies 2vi {i = x', y, z') is the same as the one 
given above at Vi, the parametric heating rate due to the 
fluctuations in the three parameters remains very small, 
Tparam < 10~*s^^. Parametric heating should thus not 
be the dominant heating mechanism in the trap. 



position of the minimum 



1 - 16e/3 4 Js 

= — — Xc- - . Zc 

3 Vl - 4£ 



1 - As 



4 _ 1 



(9a) 
(9b) 



We infer from these equations that ro is changed only 
by Xc, with a coefficient almost unity: 



Aro = 



1 



(10) 



For Xc = Q and Zc 7^ 0, the trap minima stay on the 
isomagnetic lines defined by the value rg. The cigcnaxes, 
given by the angle 9, can thus be significantly tilted for 
Zc ^ 0. The deviation A0 of the eigenaxes depends more 
on Zc than on Xr, and reads: 



4^2 _i 



/e Xc 



3(C2-l-f 3e) Vyr^^a^o 



(11) 



The dependence of the oscillation frequencies on the 
plug position is also calculated analytically at first order. 
A remarkable feature is that the highest frequency uJx' 
does not depend on the plug position at first order in Xc, 
Zc- This makes this frequency particularly stable. The 



frequency Wy 



-1/2 



depends only on Xc, as does tq. 



We express the relative shifts of ojz' and ujy as functions 
of the unperturbed position xq of a centered plug given 
by Eq. dl: 



2. Beam pointing 



Finally, we discuss the sensitivity of the trap to the 
beam pointing stability. When the plug beam is focused 
exactly in the center of the magnetic trap, the two min- 
ima are symmetrically placed with respect to the yz plane 
and present the same trap depth. However, small dis- 
placements in the plug position of the order of a few 
micrometers unbalance the depths and consequently the 
populations of the trap minima. They also modify the 
positions of the minima and the oscillation frequencies 
inside them. For large displacements, one of the minima 
can even disappear. 

Wc first remark that the symmetry with respect to 
the xz plane is still preserved with an off-centered beam, 
which implies that the trap minima always satisfy y^ = 
Moreover, the expression of the oscillation fre- 

quency ujy as a function of the gradient b' and the mag- 
netic field Bq at the trap minimum still holds. For dis- 
placements small as compared to wq, we can estimate the 
shift induced on the position of the minimum in the xz 
plane and on the oscillation frequencies. We find for the 



UJx' 

Aojj, 
Aw^ 



= 0, 



1 



1 



2 1 — e Xc 



1 Xc 

2 xq ' 



3 1 - 4e xo 3 Vl-4e xo 



(12a) 
(12b) 

-. (12c) 



From these equations, we conclude that the trap fre- 
quencies depend essentially on Xc and only slightly on 
Zc, whereas Zc essentially affects the direction of the trap 
eigenaxes. 

We can evaluate the dipolar and parametric heating 
rates from this analysis. We find again that dipolar 
excitation dominates. To limit the heating rate below 
lnK-s~^, the pointing noise of the laser at the trap fre- 
quencies should be kept below Sx = — TOdB'/im^-Hz"^ 
in the x direction, which gives a constraint three times 
stronger than in the z direction. In our experiment, a 
residual heating of 80nK-s~^ is present and is directly 
linked to measured beam pointing fluctuations of order 
Sx = — 50 dB'/xm^-Hz"^; see the Appendix. As we will 
see in the next section, this moderate heating does not 
prevent the formation of a BEC in the optically plugged 
trap. 
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III. EVAPORATION TO A BEC 



ture, the variation of T is written 



The present section is devoted to the experimental 
study of the evaporative cooling dynamics of a ®^Rb cold 
gas in the optically plugged quadrupole trap described 
in the previous section. The atoms are prepared in the 
F~l,mF = — 1 internal ground state. We first con- 
centrate on Majorana losses in a bare quadrupole t rap 
and show that the simple model introduced in Ref. 
provides a reasonable description of our experimental ob- 
servations up to a dimensionless geometrical factor that 
we are able to measure. Then, careful measurements of 
the gas phase space density during the evaporation in the 
quadrupole trap allow us to see that the cooling dynam- 
ics can be maintained in the runaway regime down to a 
temperature where Majorana losses start to prevail. We 
finally show how to optimize the effect of the plug beam 
by modifying the magnetic field gradient, to efficiently 
suppress Majorana losses and reach the BEC threshold. 



A. Measuring Majorana losses 

The Majorana loss rate r,„ of a dilute gas confined in 
a bare quadrupole trap at temperature T can be defined 
as 



It, 



h / ha \ 

= X— I I 



(13) 



where x is a dimensionless geometrical factor and k b the 
Boltzmann constant. Following the approach described 
in Ref. [l2|, one can also show that the thermodynamics 
of the trapped cloud is determined by the following set 
of equations: 



T 
T 
no 
no 



(14a) 
(14b) 



where hq is the gas peak density and Ff, the one-body 
loss rate due to collisions with the background gas. The 
temperature variation rate in Eq. (|14ap is deduced from 
the loss rate of Eq. (|T3| and from the average energy of 
atoms lost in a spin- flip event. This last average is com- 
puted under the same assumptions that led to Eq. (fT3| 
and gives an average energy of {5/2)kBT per lost atom. 
Using the energy conservation relation E = {5/2)kBTN 
and the expression of the total energy in a quadrupole 
trap E = {9/2)kBNT, it is then straightforward to ob- 
tain Eqs. ((Ti)) . 

Solving Eqs. (|14[) explicitly, one can show that no de- 
cays according to a non-exponential law involving both x 
and Ffc. In comparison, T has a simpler behavior, where 
only X enters. Denoting by To the gas initial tempera- 



where 7 = -y— 



ha 
kB 



(15a) 
(15b) 



In order to obtain the value of x experimentally, it hence 
appears more favorable to measure the gas temperature 
increase rather than the atom number decay, as it in- 
volves fewer free parameters. The heating of the gas can 
be intuitively understood as an anti-evaporation process: 
the coldest atoms of the gas are indeed the closest to the 
quadrupole magnetic field zero and thus more likely to 
undergo a spin flip. 

In order to check the validity of this model and de- 
termine the value of x for out setup, we performed the 
following experiment: After evaporatively cooling the 
atomic cloud down to T = 20 //K in a bare quadrupole 
trap with b' = 216 G-cm"^, corresponding to the maxi- 
mal value of the magnetic gradient in our setup, we de- 
compress the trap in 50 ms to a flnal value of the magnetic 
gradient 6j and keep the cloud in the trap for various 
holding times. The gas density is then measured by ab- 
sorption imaging after a time of flight of 12 ms, allowing 
us in turn to deduce T (see Sec. lIVDjl . 

The results obtained are presented in Fig. [31 and we 
can see in the inset that the gas temperature increase for 
different values of 6^ is well fitted by Eq. (|15ap . More- 
over, we correctly recover the quadratic dependence of 
the heating coefficient 7 on the magnetic gradient 6^ ex- 
pected from Eq. (jl5bp and deduce x = 0.16. This value 
agrees well with the one measured in Ref. Q for a gas of 
^•^Na atoms. This tends to prove that the simple model 
presented here describes correctly the physics of Majo- 
rana losses in a quadrupole magnetic trap, independently 
of the atomic species. 



B. Runaway evaporation in a bare quadrupole trap 

Evaporative cooling consists in truncating the highest 
part of the energy distribution of a trapped gas to force 
it to thermalize at a smaller temperature than initially. 
To describe the thermodynamics of the gas during the 
process, one usually introduces the truncation parameter 
?7, which can be defined as the ratio of the trap depth, 
fixing the maximal energy of the trapped atoms, to the 
gas temperature. It has been shown in [isj that in the 
absence of Majorana losses, the efficiency of evaporative 
cooling is completely determined by the ratio r ~ r5/Fc 
of the background loss rate F^ to the elastic two-body 
collision rate Tc- The authors of [l^ have also shown 
that below a certain value of r depending on the trap 
geometry and on 77, assumed here to be constant, the 
evaporative cooling dynamics enters the runaway regime, 
where Fc keeps increasing during the evaporation. For 
sufficiently small r, Fc and the phase-space density are 
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Figure 3: (Color online) Measurement of the Majorana heat- 
ing rate. Parameter 7 [see Eq. (|15bp ] as a function of the trap 
magnetic field gradient. The sohd hne is a quadratic fit allow- 
ing extraction of the geometrical factor x ~ 0.16, according 
to Eq. H15b|) . Inset: Cloud temperature as a function of time, 
for different magnetic field gradients, fitted by formula (|15ap 
(solid lines). See text for details. 



even expected to diverge in finite time, in turn allowing 
the gas to efficiently reach the BEC threshold. 

Majorana losses can be expected to deeply modify the 
evaporative cooling dynamics and in any case reduce its 
efficiency, especially since r„i increases with smaller tem- 
peratures. Yet maximizing Fc throughout the evapora- 
tion process still appears as the best strategy to be able 
to stay in the runaway regime during the process. This 
is achieved in our experiment by optimizing the evapora- 
tive cooling ramp shape in the bare magnetic quadrupole 
trap: the potential depth is modified through a radio- 
frequency knife whose frequency frf decreases during the 
evaporation from 50 to 4 MHz. In order to maximize Fc, 
the ramping speed v^f is increased in four steps from 3.6 
to 6.7MHz/s. Different thermodynamical quantities of 
the gas measured along the resulting evaporative cooling 
ramp are presented in Fig. |4l 

Interestingly, we see in Fig. ^a) that the gas temper- 
ature during the first part of the evaporation (region 1) 
scales as T ^ N'^-^ : this value is close to the maxi- 
mal theoretically predicted exponent of 1.2 in a linear 
trap |14| . At the same time, the phase-space density 
scales as PSD ~ N~^-^ as shown in Fig. Hljb). This re- 



sult demonstrates the efficiency of the evaporative cool- 
ing process and illustrates the weak influence of Majo- 
rana losses at these relatively high temperatures. Note 
that here the truncation parameter 77 keeps an almost 
constant value of 9. 

As soon as r < 20 fj,K, we observe a sudden decrease of 
the efficiency of the evaporative cooling, and the trapped 
gas never reaches the BEC threshold. To illustrate this 
effect we have plotted in Fig. SJ^c) the ratio Fc/Fc, which 
gives the evolution of the two-body collision rate during 
the evaporation. To compute this quantity we rely on the 
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Figure 4: (Color online) Evaporation in the bare quadrupole 
trap. Black points: experimental data for different trap 
depths (final frequency) at the maximal radial gradient of 
216G-cm~^. Red triangle: data taken immediately after 
the trap opening (see text) for a gradient of 55.4G-cm~^. 
(a) Temperature versus atom number in log-log scale. The 
blue solid line is a linear fit to the data, giving the scaling 
T ~ 7V*^'®. The dashed red line is the critical temperature 
expected for a harmonic trap with the oscillation frequencies 
given in Tabled (b) Phase-space density versus atom number 
in log-log scale, computed from the measured temperature 
and the exact knowledge of the potential shape, (c) Expected 
Fc/Fc ratio versus the atom number in linear-log scale, as 
computed from the right-hand side of Eq. H16|) . On the three 
graphs, the dashed vertical lines delimit three regions labeled 
1, 2, and 3 (see text for details). 
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simple model introduced in Ref. |13l |. which, assuming 
7/^1, predicts 

f 19 

j;^ = -Fb - yr„, + /(77)r, , (16) 

where /(ry) = | [563 + 4?7(5ry — 54)] e"''. Experimentally 
the value of Tc/Tc is estimated from the right-hand side 
of Eq. ([T6| . We see in Fig.|4l^c) that while Tc/Tc increases 
faster and faster in region 1, as expected in the runaway 
regime, it saturates in region 2 before quickly decreasing 
and reaching negative values in region 3. This study illus- 
trates how Majorana losses break down the evaporation 
efficiency and prevent the gas from reaching quantum 
degeneracy. 

C. Transfer to the optically plugged trap 

In order to circumvent the limitations imposed by Ma- 
jorana losses, it is necessary to add the optical barrier in- 
duced by the plug beam at the center of the quadrupole 
magnetic trap. This effectively decreases the atomic den- 
sity at the center of the quadrupole field, resulting in an 
exponential suppression of Majorana losses 

In our experiment, however, simply adding the plug 
beam to the bare quadrupole magnetic trap described in 
the previous section is insufficient to reduce Majorana 
losses enough to allow the gas to reach quantum degen- 
eracy. Indeed, the light shift is of order 100 ^K, result- 
ing in a potential barrier of only 50 fiK for a gradient 
of b' = 216G-cm^^, to be compared to the 20 cloud 
temperature at the point where the evaporation dynam- 
ics slows down. 

We find that it is necessary to adiabatically open the 
trap just before the gas reaches T = 20 yuK. This is done 
by reducing the magnetic gradient b' to 55.4 G-cm^^ in 
50 ms. The opening has two effects: reducing the tem- 
perature to 8 /iK and increasing the barrier to 90 /iK, 
which effectively suppresses Majorana losses. The re- 
sulting lifetime of the atoms in the opened and optically 
plugged trap reaches > 20 s. After such a sequence, r is 
reduced to ~ 10~^ which is small enough to carry on the 
evaporative cooling. 

As discussed in Sec. Ill CI the magnetic field value So 
at the trap minimum, or equivalently the trap bottom 
frequency i/q = IJ^Bq/H, may be adjusted by a controlled 
shift of the plug along the x axis, as Az/g/i^o = Xc/xq; see 
Eq. pUj) . The subsequent relative change Au)/ui in the 
mean trapping frequency ui = (waj/WyOj^/)^/^, and thus 
in the critical temperature, is about three times smaller. 
This gives room for a possible adjustment of the bottom 
frequency alone. Indeed we are able to produce quasi- 
pure BECs with trap bottom frequencies ranging from 
~ 100 kHz to ~ 300 kHz, by logarithmically sweeping z/rf 
from 2 MHz down to i/q + 50 kHz in 5 s (see the Appendix 
for details). Let us note that the smooth dependence of 
Lu with respect to Xc and Zc ensures that small long-term 



drifts in the plug position do not change much the ex- 
perimental conditions needed to reach the condensation 
threshold, apart from a shift in the final evaporation fre- 
quency. 

The results presented in this section show that the Ma- 
jorana losses can be precisely measured during the evap- 
oration process and account for the observed breakdown 
in evaporation dynamics. The temperature at which this 
occurs sets the order of magnitude of the plug barrier 
necessary to suppress Majorana losses, and an adiabatic 
trap opening might be required depending on the avail- 
able laser power. Let us note that in our setup we are 
never limited by three-body collisions, due to relatively 
low atomic densities. However, it has been shown [5| 
that trap opening also provides a way to circumvent the 
three-body losses, and thus produce large BECs in opti- 
cally plugged traps. 

IV. CHARACTERIZATION 

The hybrid optical and magnetic final trap is close to 
a linear trap at energies larger than 20 and is har- 
monic at very low energies, below 1 /iK. Between these 
two regimes, it is strongly anharmonic, and its precise 
features depend on several parameters: the magnetic field 
gradient 6', the waist and power of the plug beam, and 
its position with respect to the magnetic field zero. As 
described in Sec. [TTl the trap characteristics are mostly 
sensitive to the plug beam position. The beam position 
determines the magnetic field Bq at the bottom of the 
trap, the oscillation frequencies, and the population ra- 
tio between the two potential minima on both sides of the 
laser. In this section, we give a series of measurements 
that allow a full characterization of the trap parameters, 
and compare them to calculations of the potential. 



A. rf spectroscopy 

First, we recall that the oscillation frequency ujy, given 
by the expression ([7a|) , depends only on the magnetic gra- 
dient b' and the magnetic field at the trap bottom Bq- 
The gradient b' is known by both a calibration of the 
quadrupole coils with a Hall probe and a direct measure- 
ment with a cold cloud; see Appendix lA 11 A measure- 
ment of Bq would thus allow a prediction of the frequency 

UJy. 

The frequency vq corresponds to the rf frequency at 
which all the atoms are evaporated from the trap. Bq 
can thus be estimated directly from the evaporation pro- 
cedure. However, the rf amplitude used for evaporative 
cooling is large enough to ensure an adiabatic deforma- 
tion of the magnetic level, and thus shifts the bottom 
frequency. To get a precise value of the rf frequency at 
the bottom of the trap, we perform radio-frequency spec- 
troscopy of trapped condensates. We use an additional 
radio- frequency field linearly polarized along the y axis to 
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resonantly probe the trapped atoms. This induces spin 
flips between the F = 1, mp = —1 trapped state and the 
F ~ 1, nip ~ state in which the atoms are pushed away 
by the combination of the plug potential and gravity. We 
extract from the spectroscopy data the resonant rf field 
at the trap bottom. 

We calibrate the rf probe by monitoring the atom num- 
ber decay for small rf amplitudes at resonance. The decay 
rate is found to scale quadratically with the rf Rabi fre- 
quency, as expected from a Fermi _golden rule estimation. 
We follow the approach of Ref. [l^ to compute the de- 
cay rate as a function of the trap parameters. Although 
the purely magnetic trap used in Ref. [l^ is completely 
different from ours, the trap shape presents similarities. 
Indeed, in Ref [15|, the vertical trapping frequency is so 
small that the entire atomic cloud is displaced by the 
gravitational field to a region where the atoms experi- 
ence a magnetic field gradient. Similarly, in our setup 
the atoms are maintained in a region with a magnetic 
field gradient by the effect of the plug. In that sense, 
the magnetic landscape is the same in these two different 
traps and Eq. (7) of Ref. [l^ is still valid in our setup, 
even if the predicted spectrum width is more complicated 
to estimate, given the tilt of the eigenaxes. Thus we use 
both the measured decay rates and the measured spectral 
widths to calibrate the Rabi frequency. 

Typical spectra are obtained with about 100 Hz of rf 
Rabi frequency, as shown in Fig. [5] To avoid projection 
losses when the rf is turned on and off, we increase (de- 
crease) linearly its amplitude in 1 ms and probe the cloud 
during 100 ms. As expected Fig. [5] shows a symmetric 
spectrum with a center frequency of = 301.6±1.2kIIz. 
As seen in Sec.|TTl the oscillation frequency in the y direc- 
tion depends only on the cloud's effective distance from 
the quadrupole center tq = 27ri'o/a [see Eq. ([7a)) ]. This 
measurement of the trap bottom frequency thus gives 
To = 77.8 ± 1 iim and the resulting value of the oscil- 
lation frequency Vy = ujy/{2TT) = 76.2 ±0.8 Hz. The 
uncertainty on mainly comes from the uncertainty on 
the magnetic gradient; see Table H] 
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Figure 5: (Color online) rf spectroscopy: Number of atoms 
remaining in the trap after a 100 ms rf pulse, as a function of 
the rf frequency. The solid line is a Gaussian fit to the data. 
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Figure 6: (Color onhne) Spectroscopy of the oscillation fre- 
quencies. Squares: Measured cloud size in micrometers after 
25 ms time of flight, as a function of the modulation frequency 
of the current in the quadrupole coils. Solid red line: Multiple 
Lorentzian fit to the data. 



B. Measurement of the oscillation frequencies 

The most straightforward method of reliably measur- 
ing the oscillation frequencies of a harmonic trap consists 
in exciting dipolar and/or parametric oscillations of the 
trapped cloud. In the case of an optically plugged trap 
this can be achieved by modulating the current / in the 
coils which produce the magnetic quadrupole. In this 
way, both the cloud position and the trap frequencies are 
modulated [To| . 

In our experiment we impose I(t) = /o(l+e sin 2TTh'cxt), 
where Iq is the mean current amplitude, e = 1.8 % is the 
amplitude of the modulation and is the excitation 
frequency. After a trap modulation of 500 ms, the ex- 
citation is converted into heating during an additional 
100 ms thcrmalization time. We then switch off the con- 



finement abruptly. The atoms are finally observed after a 
25 ms time of flight (TOF) and the rms size of the cloud 
^^(^tof) is recorded. 

In Fig. |6]we display the dependence of Sz(itof) on z/ox 
measured experimentally. The different peaks in the 
graph correspond to the dipolar and parametric reso- 
nances of the trap. Their position is deduced from a 
Lorentzian fit. The first peak, centered at 76.6(4) Hz, 
directly corresponds to Uy which nicely confirms the es- 
timation of the previous section. We otherwise find 
ly^;' = 220(10) Hz and i/^/ = 121.1(5). The uncertainty 
results from the fit. Note that parametric resonances 
happen at double the frequency of dipolar resonances. 
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C. Trap parameters 

From the knowledge of the oscihation frequencies, we 
can infer the trap parameters. The magnetic gradient 
b' is measured independently to be 55.4 ± 0.6G-cm~^ 
(see the Appendix lA 1[) which means that both a and 
e = 0.0337 ± 0.0007 are known. The laser power can also 
be measured in a reliable way to be P = 5.8 ± 0.1 W by 
recording the power before and after the vacuum cham- 
ber. The three remaining independent parameters, wq, 
Xc, and Zc, which are more difficult to obtain by a di- 
rect measurement, are deduced from the three measured 
oscillation frequencies. 

The value of ujy directly gives tq, which is the isomag- 
netic surface on which the minimum lies, from Eq. (|7ap. 
which also holds for an off-centered plug beam. The same 
information also comes from the rf spectroscopy. Now, 
the general idea for finding the trap parameters is to first 
adjust the tentative value of wq to fit the correct lJx' ■ As 
does not depend on the plug position at first order, 
the expression for ^ Q, which enters in ujx' , can be safely 
used. Let us introduce the light shift uq = Uo/{hawo) 
in dimensionless units, uq is a function of P and wq- 
Mo = kP/wQ, where fc is a fixed parameter for a given 
gradient and is equal to fc = 1.98 x 10^ • /im"^ in our 
case with the gradient of Table HI The function ^(wq) is 
thus known for the laser power P = 5.8 W, which allows 
us to determine wq from ujx' by inverting the relation 



(17) 



We find Wo = 46 /im, in fair agreement with an estimate 
from a direct optical measurement. 

The knowledge of wq allows prediction of a zero-order 
value Tq for the effective radius from Eq. From 
the shift Aro = ^o — Tq of this value compared with 
the measured one, we deduce the horizontal plug shift 
Xc = Vl — 4eAro; see Eq. ([T0| . Finally the value of Zc 
is chosen to fit w^- , a first guess being given by Eq. (|12cp 
where the zero-order frequency comes from Eq. (j7cp . Ap- 
plication of this method allows us to determine the values 
of Wo, Xc, and Zc given in Table U] and calculate the trap 
depth. The uncertainties are deduced from the relations 
(fTT]) . ([5|), and (fT2|) and from the uncertainties of the ex- 
perimental measurements. 



D. Time-of-flight analysis 

In order to access the temperature of the trapped 
atomic clouds wc rapidly switch off the confinement 
(< 0.5 ms) and let the gas expand during a controllable 
duration <tof- Analysis of the shape of the atomic col- 
umn density obtained from absorption imaging (see the 
typical data displayed in Fig. [7]) allows us to deduce the 
temperature of the trapped gas. Depending on the phys- 
ical regime of the cloud, different methods are used. 



A simple strategy to obtain the temperature of the gas 
consists in measuring the rms width of the cloud tT„ , with 
u = x,y, z. For a freely expanding gas it can be expressed 
as 



(itof ) = al (0) + 



Hoi 



where 



al (0) = C y dr exp[-/3;7 (r)] 



(18) 



(19) 



with /3 = \kBTY^ and C = [/ dr exp[-/3;7 (r)]] . For 
sufficiently high temperatures, the effect of the optical 
plug on the trapped gas can be neglected and U (r) re- 
placed by J7b (r) -I- Mgz in Eq. ([T^. The calculation of 
a\ (0) then becomes analytical and using Eq. it is 
straightforward to deduce T from a single picture. As 
soon as the gas gets too cold, approximating U (r) by 
Ub (r) + Mgz is no longer valid and a\ (0) can be ob- 
tained only numerically. At this point, a better strategy 
to deduce T consists in taking a series of images with 
different itof and use Eq. (fT5|) . 

Below the Bose-Einstein condensation threshold, the 
expansion of the gas cannot be considered as ballistic and 
we have to rely on a different strategy. We hence fit an 
analytical formula based on ideal Bose gas theory [l^ to 
the thermal tail of the atomic density distribution and de- 
duce the value of T. In principle, such an approach could 
be extended to the investigation of the BEG parameters, 
like the initial Thomas- Fermi radii [T^j. However, tech- 
nical limitations of the imaging system have prevented 
us from applying this analysis. 



V. CONCLUSION 

In this paper, we describe the successful evaporation of 
a rubidium 87 gas to Bose-Einstein condensation in a lin- 
ear quadrupole trap plugged by a laser beam at 532 nm. 
We provide a simple model which describes the trap char- 
acteristics. We check its predictions quantitatively by a 
direct measurement of the oscillation frequencies. For ex- 
ample, the value of Uz' /(^Jy = 1.61 predicted by the sim- 
plified model with a centered plug [see Eq. (fZc)) ] agrees 
with a direct measurement within 2%. 

We present two spectroscopic tools to fully character- 
ize the trap. Other possible diagnostics include the sep- 
aration between the two minima, measured by in situ 
absorption imaging with the laser and the magnetic field 
on (see the inset of Fig. ^ or the orientation of the 
eigenaxes measured after a free flight. These methods 
could improve the error bars on the measured param- 
eters. However, the two proposed methods measure a 
frequency directly and are thus more reliable as they do 
not require a particularly high imaging quality. In partic- 
ular, the measurement of the bottom frequency gives 
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Figure 7: (Color online) Analysis ol time-of-flight image, (a) 
Absorption image of a cold atomic cloud after 25 ms time of 
flight. High optical density at the center of the image indicates 
the presence of a BEC. Thermal tails are yet still visible, (b) 
Integrated profile of (a) along the x axis. The dashed blue 
line corresponds to the integrated profile of a two-dimensional 
(2D) fit of the thermal tail of the atomic cloud displayed in 
(a), (c) Same as (b) with integration along the z axis. 



access to the main scaling frequency LOy. Provided an in- 
dependent estimation of the plug beam waist is available, 
this single measurement allows a good estimation of the 
average frequency w, which is the natural scale for the 
critical temperature or the chemical potential [l8l |. 

A detailed study of the influence of the plug position 
on the trap parameters confirms that the trap character- 
istics are robust with respect to small deviations in the 
experimental parameters. The only relevant source of 
heating comes from beam pointing fluctuations. Again, 
their effects remain modest. Without active stabiliza- 
tion of the pointing, the dipolar heating rate is kept to a 
reasonable value of 80nK-s~^ in our experiment; see the 
Appendix. 

Evaporation in the optically plugged quadrupole trap 
is shown to lead efficiently to degeneracy, with a dy- 
namics comparable to evaporative cooling of rubidium 
in similar linear traps Q. Another result of the paper is 
the measurement of the effective volume entering in the 
model of Majorana losses [ll|. We find a value in agree- 
ment with recent measurements in a similar setup with 
sodium atoms Q. 

Finally, we have shown that a quadrupole trap op- 
tically plugged with a 532 nm focused beam, initially 
demonstrated with sodium atoms [J] , is also well adapted 
to the production of degenerate rubidium gases, despite 
the much larger detuning. It should also work as well 
with atoms with intermediate wavelengths for the main 



dipole transition, like lithium (671 nm) or potassium 
(767 nm). This makes this trap particularly well suited 
to production of mixtures of degenerate gases. 
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Appendix A: Experimental setup 

In this appendix we first give an overview of the differ- 
ent elements of the experimental setup. Then we describe 
the laser sources and finally give technical details on the 
experimental sequence. 



1. Overview of the experimental setup 

The experimental setup can be decomposed into three 
main parts as depicted in Fig. [5] A first vacuum cham- 
ber is filled with a hot vapor of ^""Rb atoms heated in 
an oven to 70 °C. This sets the ®''Rb partial pressure to 
10~^mbar. The atoms are collected by a bi-dimensional 
magneto-optical trap (2D MOT) [l^. An additional laser 
beam pushes the atoms into a second vacuum chamber 
where they are captured by a three-dimensional magneto- 
optical trap (3D MOT). We typically load about 2.5 x 10^ 
atoms in the 3D MOT in 5 s, these numbers being mostly 
sensitive to the power and direction of the pushing laser 
beam. Thanks to a differential pumping tube the partial 
pressure in the 3D MOT chamber can be kept of the or- 
der of 10^^° mbar. as deduced from the 3D MOT lifetime 
of 30s [23. 

One of the main features of our setup comes from the 
water-cooled 3D MOT quadrupole magnetic field coils 
which are held on a motorized translation stage that can 



be displaced at will along the x axis [21| . After the trans- 



fer of the atoms from the 3D MOT into a quadrupole 
magnetic trap made with the same coils, the translation 
stage is moved along a 28.5 cm path, bringing in turn a 
fraction of the atoms into a 10 x 10 mm^ inner-size science 
glass cell (Starna) with very good optical access. The 
cell walls have a width of 1.25 mm and are anti-reflection 
coated on their external sides at 532 nm and between 650 
and llOOnm. During the transport in the moving trap, 
the atoms pass through a 4-mm-diameter, 94-mm-long 
tube, which ensures an ultrahigh-vacuum environment 
in the flnal chamber. The atoms are flnally transferred 
into a second quadrupole magnetic trap induced by a pair 
of water-cooled conical coils made of 40 loops of hollow 
copper tubes each. There the loss rate Fb due to colli- 
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sions with the background gas corresponds to a hfetime 
= 120s. 

The cahbration of the magnetic gradient is done di- 
rectly with an ultracold cloud, in the foUowing way. The 
magnetic field is suddenly switched off for a short du- 
ration, such that the atomic cloud starts to fall in the 
gravitational field. It is then switched on again, and the 
atoms oscillate along the z vertical fiber axis in both the 
quadrupole and gravitational fields. The two gradients 
add when the cloud is above the magnetic zero and sub- 
tract below. From a parabolic fit to the data and the 
knowledge of <?, we infer a calibration of the imaging sys- 
tem and the value of the magnetic gradient. 



2. Laser system 

All the 780 nm light sent onto the atoms is prepared on 
a separate optical table and carried through single-mode 
polarization-maintaining fibers. 

For cooling and pushing Rb atoms in the 2D and 3D 
MOTs, we built an agile and powerful 780 nm laser source 
relying on frequency doubling of an amplified Telecom 



laser [22|. A fiber-pigtailed, 2mW, distributed-feedback 



laser diode (Fitel) emitting at 1560 nm feeds a 40-dB er- 
bium doped fiber amplifier (Keopsys), with a maximum 
output power of 10 W. The second harmonic is gener- 
ated in a 3 X 0.5 x 50 mm"^ periodically poled lithium nio- 
bate (PPLN) crystal (HC Photonics). The quasi-phase 
matching condition is met by regulating the tempera- 
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Figure 8: (Color online) Scheme of the experimental setup. 
The three main chambers (2D MOT, octagonal chamber for 
the 3D MOT and science cell) are visible, as well as the rail for 
the transport of the magnetic coils (on the left). The MOT 
beams propagate at 45° with respect to the x and y axes. The 
moving coils are represented by a blue hollow disk. The plug 
beam (in green) is aligned with the y axis. The x and y axes 
cross at the center of the quadrupole trap. The remaining 
part of the vacuum chamber (ionic pumps, etc.) have been 
omitted for clarity. 



ture of the crystal at 85 ± 0.1 °C with a home-made 
oven. The 1560 nm beam has a 70 //m waist to maximize 
the doubling efficiency according to the Boyd-Kleinman 
model [2^. In single pass, we obtain a maximum of 2 W 
of linearly polarized light at 780 nm in a TEMOO mode. 

The frequency control of the doubled laser is made by 
beat note locking with a reference laser. This reference 
laser is a 780 nm, 70 mW, narrow linewidth external cav- 
ity laser (RadiantDye) locked by the saturated absorp- 
tion technique. Its free linewidth has been measured to 
3 kHz and its linewidth when locked has been estimated 
at around 180 kHz. The beat note between the reference 
laser and the doubled laser is frequency locked to an rf 
signal of adjustable frequency around 270 MHz. 

By tuning the beat-note frequency we are able to sweep 
the doubled laser frequency over a large span from -1-1 OF 
to -60F around the 551/2,^ = 2 ^ ^P?,/2,F' = 3 cy- 
cling transition of the ®''Rb D2 line, where F is the tran- 
sition linewidth, without altering the output intensity. 
The intensity is independently controlled or switched ofl[ 
by an acousto-optic modulator. The doubled light is then 
split into four beams injected into polarization maintain- 
ing fibers. Two transverse cooling beams and a weak 
pushing beam, with a total power of 120 mW, seed the 
2D MOT. The last one is used for the 3D MOT and is 
coupled to one of the two input ports of a 2 x 6 fiber 
cluster (Schafter-|-Kirchhoff). 

At the cluster output, each of the six fibers is connected 
to a compact three-lense collimator system (SYRTE Labs 
design). The coUimated beams are clipped to a diame- 
ter of 1 in. by a quarter-wave plate directly set at the 
collimator output to produce the required circularly po- 
larized light. The total intensity of the six cooling beams 
is 41 mW-cm^^. 

The repumper laser is a 70 mW Sanyo laser diode, fre- 
quency locked to the F ^ 1 ^ F' — 2 transition of the 
D2 line. It is superimposed onto the 2D MOT transverse 
beams before their injection into the fibers, while it is 
mixed with the 3D MOT beams through the second in- 
put of the fiber cluster. The repumper laser is not present 
in the pushing beam, which limits the 2D MOT atomic 
beam velocity and improves the recapture efficiency [l^] . 

The plug beam originates from a high-power diode- 
pumped laser at 532 nm (Spectra Physics Millennia X) 
with an output power reaching 10 W in a TEMqo spatial 
mode. To control the laser intensity and allow for the fast 
switching of the beam, we use an acoustic-optic modula- 
tor with a center frequency of 110 MHz. At the position 
of the atoms, the optical plug power is about 6W for a 
waist of 46 /im. The beam position is controlled by two 
piezoelectric actuators on a mirror mount. 

The power spectrum density of intensity fluctuations of 
the laser beam amounts to — 93dB-Hz^^ at the trap fre- 
quencies, which corresponds to a calculated heating rate 
of about InK-s"^; see Sec. Ill Cl Using a quadrant pho- 
todiodc, we characterized the pointing stability of the 
plug beam. The long-term beam pointing stability is 
rather good, with a drift below 1 /im over 1 week. We 
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also recorded the power spectrum of beam pointing fluc- 
tuations. We find < — 50 dB-^m^-Hz~^ over the trap 
frequency range. This value agrees with a measured heat- 
ing rate of 80nK-s~^ in the optically plugged quadrupole 
trap. This figure could be improved if necessary by ac- 
tively locking the plug position to a reference measured 
with the quadrant photodiode. Finally, the waist fluctu- 
ations typically take place at very low frequencies, linked 
to thermal effects in the laser cavity. We did not quan- 
tify this level, but we expect an extremely low value at 
the trap frequencies, with negligible contribution to the 
heating rate. 



3. Experimental sequence 

Before the transfer of the atoms from the 3D MOT 
to the quadrupole magnetic trap, the 3D MOT is 
first compressed by progressively increasing the mag- 
netic field gradient from 9.5 to 45G-cm~^ and increas- 
ing the 3D MOT laser beams detuning from — 3r to 
— 14r. The laser beams are then turned off and the mag- 
netic field axial gradient is linearly increased from 45 to 
200G-cm~^ in 250 ms, allowing the capture in the mag- 
netic field quadrupole trap of about 10^ atoms in the 
\F = l,mF = —1) state. 

The displacement of the quadrupole trap from the 3D 
MOT chamber to the science glass cell relies on a motor- 
ized translation stage (Parker 404- XR) whose maximal 
velocity, acceleration, and jerk are controlled by built- 
in software. During the first half of the transportation 
the maximal acceleration of the translation stage is set 
to 0.8m-s~^ and its jerk to SOm-s"'^. It hence reaches 
a maximal velocity of Im.s"^. The second half of the 
transfer corresponds to the mirrored image of the first 
half. Finally the transfer efficiency is about 20% and 
is mostly limited by the free evaporation of the hottest 
atoms against the walls of the small tube during the dis- 
placement of the translation stage. 

As soon as the atoms have reached the science glass 
cell, the current in the transfer magnetic field coils is 
ramped down to zero while the current in the conical 
magnetic field coils is ramped up from zero in order to 
keep the magnetic field gradients constant. Up to 3 x 10^ 
atoms with a temperature of 150 fxK are loaded in the 
final magnetic field quadrupole trap. The current is pro- 
vided by a Delta Elektronika power supply, which has a 
low relative current noise, — 125dB-Hz~^ at most, and 
the associated dipolar heating is negligible. In order to 
increase the collision rate, the axial magnetic field gra- 
dient is finally adiabatically ramped up to 432G-cm~^, 
the radial gradient thus being b' = 216G-cm~^. In this 
trap the atom lifetime is larger than 120 s at the initial 
temperature of 270 ^K. 

The plug is then switched on and the rf-induced evap- 
oration starts. The plug alignment proceeds as follows. 
We first leave the plug on during the time of fiight and 
absorption imaging which allows us to image the plug 



as a depression in the expanding gas density profile. In 
this way, we are able to roughly match the plug posi- 
tion to the magnetic trap center (known from in situ 
absorption imaging). Then, at lower temperatures (typ- 
ically at 20 /iK), the plug position is optimized to obtain 
a symmetric expansion of the cloud due to initial acceler- 
ation induced by the plug during the time of flight. Once 
this optimization is done, we proceed with the evapora- 
tion sequence and the adiabatic trap opening discussed 
in Sec. IIIII Finally we switch off the plug and the mag- 
netic trap simultaneously and optimize the plug position 
on the peak density after a given time of flight. We re- 
peat this for lower and lower rf frequencies, until the 
expanded cloud density p rofile starts to deviates from a 
pure Gaussian. At this point BEG is observed on further 
lowering of the RF frequency. Eventually we finely tunc 
the plug position to place the BEG on a given magnetic 
isopotential. 

Both the rf field applied for evaporation and the prob- 
ing rf field are linearly polarized along the y axis. The 
evaporation field is generated by the rear panel out- 
put of a direct digital synthesizer (DDS) (TaborElec 
WW1072), going through a 5 W amplifier (Kalmus) and 
roughly impedance matched to a 10-mm-radius 20-loop 
coil, to avoid any unwanted resonance in this frequency 
range. The overall circuit length, including rf switches, 
is much smaller than the shortest wavelength to avoid 
high-frequency resonances. 

The frequency is ramped from 50 down to 4 MHz with 
a piece-wise linear time dependence, in order to keep ef- 
ficient evaporation dynamics in 13.5 s. At this point we 
have about 10^ atoms at 20 /iK. The trap is then adi- 
abatically opened to a 55.4G-cm~^ gradient in 50 ms, 
which results in a cloud temperature of 8 /jK. After the 
trap opening, the evaporation is carried on with a second 
antenna (30 mm radius, 20 loop coil, about 30 mm away 
from the trap center), producing a rf field linearly po- 
larized along the y axis. This antenna is fed directly by 
a Stanford DS345 function generator, producing a loga- 
rithmic ramp from 2 MHz to 300 kHz in 5 s. The probe 
antenna for radio-frequency probing the condensate is a 
30-mm-radius 20-loop coil, directly fed by another Stan- 
ford DS345. 

We use absorption imaging to measure the density dis- 
tribution of the atoms once released from the optically 
plugged quadrupole trap. The probe beam is derived di- 
rectly from the reference laser, its waist at the atoms 
position is 1.6 mm for a total power of 300 /zW. It is 
circularly polarized and is superimposed onto the plug 
beam along the y axis thanks to a dichroic mirror. Dur- 
ing the 100 fis imaging pulse, a bias magnetic field of 
1 G is applied along y to define a quantization axis. Af- 
ter transmission through the science cell, the probe and 
plug beams are separated with a polarizing beamsplitter 
(at 532 nm, transparent and polarization independent for 
the probe beam at 780 nm). The cloud is imaged with 
a single lens of focal length 250 mm onto a GGD cam- 
era (Andor IXON-885D) with a pixel size of 8 fim. The 



13 



magnification is 0.94, as measured by the procedure de- came ra pixels by the remaining optical plug photons, we 
scribed in Appendix lA II To avoid the saturation of the use a green- light filter (Layertec). 
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